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Two-level modeling for nanoscale pattern formation on silicon target by Ar+ ion sputtering is presented. Phase 
diagram illustrating possible nanosize surface patterns is discussed. Scaling characteristics for the structure 
wavelength dependence versus incoming ion energy are defined. Growth and roughness exponents in differ- 
ent domains of the phase diagram are obtained. 
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1. Introduction 

It is well known that low and medium energy ion sputtering may induce a fabrication of 
periodic nanoscale structures on an irradiated surface [1]. Depending on the sputtered substrate 
characteristics and sputtering conditions, different types of nanoscale structures such as ripples, 
nanoholes and nanodots can grow on a target during ion beam sputtering [2-6]. These patterns have 
been found on both amorphous and crystalline materials including insulators, semiconductors and 
metals (see reference [7| and citations therein). Main theoretical models describing ripple formation 
are based on the results of famous works by Bradley and Harper [4], Kardar et al. [8], Wolf and 
Villian [9], and Kuramoto et al. [10]. The main control parameters in these models reduced to 
surface tensions, tilt-dependent erosion rates and diffusion constants are determined by sputtered 
substrate characteristics and sputtering conditions (see for example [7]). 

Among theoretical investigations there are a lot of experimental data manifesting a large class 
of patterns formed due to the self-organization process. It was experimentally shown that the main 
properties of pattern formation processes depend on ion-beam parameters such as ion flux, energy 
of deposition, angle of incidence, and temperature of the substrate (target). Therefore, to study 
the ion beam sputtering processes theoretically one needs to determine the mentioned parameters 
of the model according to the physical conditions related to concrete materials. 

One of the most frequently used materials for ion beam sputtering is silicon because it is the 
mainstream material in modern microelectronic industry and it is readily available with high purity 
and quality. Nanostructuring of silicon has received much attention due to its potential application 



in developing the Si light sources 11]. Various techniques such as acid etching, ion implantation, 
reactive evaporation, chemical vapor deposition and molecular beam epitaxy have been used in 
developing the Si nanomaterials (porous Si, Si nanocrystal-doped dielectrics and Si quantum dots) 
(see reference [ll| and citation therein). 

In this paper we study the properties of the formation of nanoscale patterns on a silicon target 
sputtered by Ar + ions. To this end, we use a two-level scheme, based on Monte-Carlo simulations 
and the modified Bradley-Harper theory. In the first approach we compute the ion energy dependent 
penetration depth, widths of the ion energy distribution and sputtering yield. Next, we exploit these 
characteristics as input data for the continuum approach describing the evolution of the surface 
height field. We define the domains of values for the angle of incidence and ion energy where 
different nanoscale structures can be formed. The dynamics of nanoscale pattern formation is 
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discussed. We will show that at fixed values for the incidence angle, one has two scaling exponents 
for wavelength related to small and large values for ion energy. In addition, we obtain roughness 
and growth exponents. 

The work is organized in the following manner. In section 2 we present the theoretical model in 
the framework of the modified Bradley-Harper approach. In section 3 using Monte-Carlo modeling 
we compute the main characteristics incorporated into the continuum theory. The related phase 
diagram, the dynamics of the formation of nanoscale structures, the ion energy dependent wave- 
length and scaling properties of the surface morphology are discussed in section 4. We conclude in 
the last section. 



2. Theoretical model 

Let us consider a d-dimensional substrate and denote with r the d-dimensional vector locating 
a point on it. The surface is described at each time t by the height z = h(r,t). If we assume that 
the surface morphology changes during ion sputtering, then we can use the model for the surface 
growth proposed by Bradley and Harper [4] and further developed by Cuerno and Barabasi [2]. 
We consider the system where the direction of the ion beam lies in x — z plane at an angle 6 to the 
normal of the uneroded surface. Following the standard approach one assumes that an averaged 
energy deposited at the surface (let say at point O), due to the ion arriving at the point P in the 
solid, follows the Gaussian distribution Q E(r) = [e/(27r) 3 ^ 2 ap 2 ] exp [—z 2 /2a 2 — (x 2 + y 2 )/2p 2 ] ; 
e denotes the kinetic energy of the arriving ion, a and \i are the widths of distribution in directions 
parallel and perpendicular to the incoming beam. Parameters a and \i depend on the target material 
and can vary with physical properties of the target and incident energy. The erosion velocity at 
the surface point O is described by the formula v = p f n dr<£(r)£?(r), where integration is provided 
over the range of the energy distribution of all ions; here 3>(r) are the corrections for the local slope 
dependence of the uniform flux J. The material constant p is defined as: p = 3/ (^^(NUqCq) -1 , 
where Uq and Co are the surface binding energy and the constant proportional to the square of the 
effective radius of the interatomic interaction potential, respectively (l2|. The general expression 
for the local flux for surfaces with non-zero local curvature is [13]: 



y, h) = J cos 



arctan 



^(v x h) 2 + (v y hy 



Hence, the dynamics of the surface height is defined by the relation d t h ~ —v(Q — V x h, V 2 /i, V 2 /i) 
and is given by the equation d t h ~ -v(0)y/l + (V/i) 2 , where < < tt/2 0-0, H, [3]. The linear 
term expansion yields d t h = — + J^xh + ^V^/i; where V = <9/<9r, V a = d/da, a = {x,y}. 
Here vo is the surface erosion velocity; 7 = 7(0) is a constant that describes the slope depending 
erosion; v a = u a (9) is the effective surface tension generated by erosion process in a direction. 

If one assumes that the surface current is driven by differences in chemical potential /i, then 
the evolution equation for the field h should take into account the term — V • j s in the right hand 
side, where j s = KV(V 2 h) is the surface current; K > is the temperature dependent surface 
diffusion constant. If the surface diffusion is thermally activated, then we have K = D s K P /n 2 T, 
where D s = Do exp (—E a /T) is the surface self-diffusivity (E a is the activation energy for surface 
diffusion), n is the surface free energy, p is the areal density of diffusing atoms, n is the number 
of atoms per unit volume in the amorphous solid. This term in the dynamical equation for h is 
relevant in high temperature limit which will be studied below. 

Assuming that the surface varies smoothly, we neglect spatial derivatives of the height h of 
third and higher orders in the slope expansion. Taking into account nonlinear terms in the slope 
expansion of the surface height dynamics, we arrive at the equation for the quantity h! = h + v$t 
of the form 0, 0] 

d t h = jV x h + VcSlji + X - \ a (V a h) 2 - KV 2 (V 2 h) + f (a, 2/, t), a = {x, y}, (1) 

where we drop the prime for convenience. Here we introduce the uncorrelated white Gaussian noise 
£ with zero mean mimicking the randomness resulting from the stochastic nature of the ion arrival 
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to the surface. In equation (pQ) the effective surface tensions v x and v y generated by the IBS, the 
tilt-dependent erosion rates \ x and \ y are defined through the incident angle penetration depth 
of incident ion a, distribution widths a, fi and the sputtering yield Yq as follows (2,01: 

7 = F j^ [a\a\(? (a* - l) - a^s 2 ] , 

Oj (j / A A 4.999 9999 4 4\ 

v x = Foa-^ (2a a s - a a a^s c + a a a^s c - a^c ) , 



-F a 



c 2 a 2 



2/ ' 



Fo^j- 4 Kay (3 + 2c 2 ) + Malays 2 - a^c 4 (l + 2 S 2 
- f 2 [2aU 2 ~ ay (l + 2 S 2 )] - a^cV - / 4 } , 



X y = Fo^iay+alay-aUy-f). (2) 



\± V 7J o n 

4 2 i „2„2 2 „4 2 2 ^2\ 



2/ 

Here we have used the following notations: 



F o = 7=^ exp —f— , (3) 




f = a a s +a c . (4) 

Let us perform the stability analysis for a system with additive fluctuations. To this end, we 
average the Langevin equation (pQ) over noise and obtain 

dt (h) = 7 V a (h) + v x V 2 xx (h) + v y V 2 yy (h) + ^ ((V^) 2 } + ^ ((V^) 2 } - ifV 4 (h) . (5) 

Considering the stability of the smooth surface characterized by (h) = 0, we can rewrite the 
linearized evolution equation in the standard form: 

d t (h)= (f>ef + £ e f) (h), (6) 

with notations 

i>ef = 7V X + v x V 2 xx + v y W 2 yy , K ei = -KV 4 . (7) 

It is easy to see that equation ([7]) admits a solution of the form (h) = A exp[i(k x x+k y y—ujt)+xt\ • 
Indeed, substituting it into equation ([7]) and separating real and imaginary parts we get 

w = -l(0)k x , (8) 
X = -v x {9)k 2 x -v y {e)k 2 y -K{k 2 x + k 2 y ) 2 . (9) 

As far as Fq , /, a are positive values, hence one has v y < 0, whereas v x can change its sign. 
Therefore, the Bradley-Harper model does not provide for stable smooth surface. Hence, we can 
conclude that if v x > 0, then ripples (wave patterns) appear in x-direction. On the contrary, when 
v x < 0, equiaxed structures (nanodots/nanoholes) can be formed on an eroded surface. In addition, 
the sign of the product • X y can play a crucial role in ripple formation processes |15| . 

For the noiseless nonlinear model (pQ) it was shown that as the sets v a and X a are the functions 
of the angle of incidence 6 G [0,7r/2] there are three domains in the phase diagram (a a , 0) where 
v x and \ x change their signs, separately [2]. This results in the formation of ripples in different 
directions x or y varying a a or 0. 

One needs to note that in the Bradley-Harper approach describing the processes of ripple 
formation on amorphous substrates the penetration depth can be approximated as a(e) ~ e leading 
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to the power-law asymptotics for the wavelength of the ripples A versus ion energy as follows: 
A ~ e -1 / 2 [4]. In the next section, performing calculations for the sputtering of the silicon target 
by Ar + ions we shall show that there are deviations from these asymptotics due to power-law 
dependence of a, a, \i versus ion energy e. Moreover, it will be shown that the sputtering yield 
depends on both incident angle end ion energy e in a power-law form. To define a, a, /i and the 
sputtering yield, we shall use Monte-Carlo approach. 



To study the evolution of the silicon surface morphology during Ar + ion beam sputtering one 
needs to know the energy of ions and target characteristics such as: penetration depth of the Ar + 
ions into the silicon target a; widths of the distribution in parallel and perpendicular directions of 
the incoming beam (<r and /i) and sputtering yield Yq . Moreover, to simulate the target morphology 
evolution one should define temperature T, uniform flux J, atomic density of the target N, surface 
binding energy Uo and the effective radius of interatomic interaction potential. Using data from 
reference [5] at T = 550 C one has K = C'8.49 x 10 3 nm 4 /s, where addimer concentration 
is C = 0.04 atoms/site (~ 4% coverage), or, C = 0.07 atoms/nm 2 . To define the material 
constant p we shall use TV ~ 50 atoms/nm 3 , Uo = 4.73 eV. For the effective radius of interatomic 
interaction potential, we put the length of the main diagonal of silicon primitive cell with lattice 
parameter 0.5437 nm. To compute the time evolution of the silicon surface morphology, we should 
calculate the effective surface tensions v x and v y generated by the IBS, and the rates and X y 
[see equation @]. Following relations ([3j) and dependent on ion energy, parameters a, a and //, 
as far as Yq = Yq(0, e), can be computed. From experimental point of view, the control parameters 
at IBS are the energy of ion beam, off-normal incidence angle and ion flux. In all our calculations 
we put J = 20 ions/(nm 2 s). Thereafter we vary the ion energy in the interval 100 eV-i-10 keV and 
use intermediate off- normal incidence angles 6 G [40°, 65°]. 

In further study, we use the well-known program codes (TRIM and SRIM) to calculate the 
stopping range of ions in matter and transport range of ions in matter. Description of algorithms 
and the basic principles for Monte- Carlo calculation of both the transport range of ions in matter 
and the stopping range of ions in matter can be found in jTgj; TRIM and SRIM codes can be found 
on the web-site |www.srim.org| 

Values for parameters a, a and \i for silicon target sputtered by Ar + ions were obtained with 
the help of SRIM code (a program for calculating the stopping range of ions in matter). Results 
for relative penetration depths a G = a/ a and a M = a/p versus ion energy are shown in figure [T] 
(dependencies a(e), a(e) and n(s) are shown in the insert). In figure [T] it is seen that longitudinal 
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Figure 1. Relative penetration depths of the Ar + ions into the silicon target versus ion energy 
(penetration depth and distribution widths versus ion energy are shown in insertion). 
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and transverse widths a and /i, respectively, as far as a a and a M satisfy the following relations, are 
as follows: a < \i and a a < a M . In reference [17| it was shown that the rotated ripple structures 
formed when A^A^ < with rotation angle cp = tan -1 \/—X x /X y can be observed at small incidence 
angles when a M < a a (a a = 1) and at intermediate and large when a M > a a . Hence, one can 
expect the appearance of rotated ripple structures in our system. 

It is principally important that dependencies of penetration depth a and longitudinal and 
transverse widths a and /i, respectively, versus ion energy e deviate from the linear law predicted 
by Bradley and Harper For a silicon target sputtered by Ar + ions we have obtained a power-law 
approximation of the form: <p(e) = A-\- Be c , where <fi = {a, a, /i}, constants A : B and C are fitting 
parameters. So, we can expect that the wavelength dependence A(e) ~ e~ 6 can be characterized 
by the exponent 5 ^ 1/2. In our further continuum approach we shall use the obtained power-law 
asymptotics for a, a a and a M from Monte- Carlo simulations. 

To compute the dependence of the sputtering yield versus ion energy and angle of incidence we 
use Monte- Carlo approach realized in TRIM code (program for the calculation of transport range 
of ions in matter). The results of calculations for sputtering yield versus incident angle at fixed 
ion energy and sputtering yield versus ion energy at fixed incident angle are shown in figures [2] (a) 
and[2](b), respectively. In figure [2] it is seen that sputtering yield depends on both ion energy and 
incidence angle in accordance with a power law as follows YbWO = A! + B f ip c , where i/j = {0,s}, 
constants A' B' and C are fitting parameters. Therefore, all parameters (y x , v y , X x , X y ) required 
to monitor the time evolution of silicon surface morphology during IBS are well defined. 




4. Surface morphology change during sputtering 

4.1. Phase diagram and typical patterns 

Firstly, let us compute a phase diagram e{6) defining domains for different surface patterns of 
silicon sputtered by Ar + ions. To this end, we shall monitor a sign change of surface tension v x and 
tilt-dependent erosion rates X x and X y (as it was mentioned above v y is always less than 0). The 
corresponding phase diagram indicating possible patterns is shown in figure [3] We need to stress 
that in the related interval for both the ion energy and the incidence angle except v y < 0, one has 
X y < 0. From figure [3] it follows that plane (0,e) is divided by three curves into five domains A, 
B, C, D and E. If one crosses the dash-dot curve, then quantity v x changes it sign. Therefore, in 
the linear regime at small incidence angles 6 (domain A), instability of the silicon surface occurs 
in both x and y directions due to v y < and v x < 0. In the domain E (at large 6) in the linear 
regime, patterns are stable in ^-direction due to v x > 0. At large times (nonlinear regime) the 
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surface morphology is governed by nonlinear parameters A^ and X y . Solid curve in figure [3] divides 
domains characterized by A^ < and X x > 0. Therefore, between solid and dash-dot lines only 
X x is positive (domains C and D), whereas in the domain E both v x and A^ are positive. Dash 
curve corresponds to the condition v x = v y . Hence, before the dash curve (domains A and C) 
when v x < v y: vertical elongated surface structures should be formed, whereas after the dash curve 
(domains B, D and E), the corresponding structures should be of a horizontal elongated type. 




(cleg) 



Figure 3. Phase diagram and typical surface patterns. 

To illustrate typical structures in each domain in figure [3] we numerically solve equation (pQ) 
on quadratic lattice L x L of the linear size L = 256 with periodic boundary conditions. Spa- 
tial derivatives of the second and fourth orders were computed according to the standard finite- 
difference scheme; the nonlinear term (Vh) 2 was computed according to the scheme proposed 
in references [HI, ED]. We have used Gaussian initial conditions taking (h(r,t = 0)) = and 
((Sh) 2 ) =0.1; the integration time step is At = 0.005 and the space step is t = 1. 

Typical surface patterns in domains (A-E) are shown in figure [3j It is seen that on the left hand 
side of the solid curve when v y < and v x < 0, pattern type of holes is realized (see snapshots 
A and B). It follows that patterns realized at high energy ions are characterized by small size 
(see snapshot A), whereas at small e one has large-scale patterns (see snapshot B)Q|. Moreover, 
orientation of holes in points A and B is different. It is defined by a minimal value of both v x 
and v y . Structures, shown by snapshots C and D (ripples) are characterized by positive value of 
parameter A^ , which defines nonlinear effects in x-direction. An orientation of the corresponding 
ripples is defined by a minimal value of both v x and v y as in the previous case. Hence, as far as 
v x < v y from the left of the dashed curve in figure [3j the related patterns in domains A and C are 
elongated in y direction. On the contrary, in snapshots B, D, and E there are horizontal elongated 
structures. Structures in snapshots A, B, C and D are characterized by instabilities in both x and y 
directions due to v x < and v y < 0. In the domain, indicated by point E due to v x > 0, structures 
are stable in x direction. 

The obtained phase diagram is in good correspondence with the results of experimental studies 
of the dynamics of the surface Si(001) sputtered by Ar + ions [20|, where according to the exper- 
imentally obtained a phase diagram in the plane "ion energy - angle of incidence" it was shown 
that if the angle of incidence or ion energy varies, then orientation of ripples can be changed. 

1 Dependence of wavelength versus ion energy at fixed values for incidence angle will be discussed later. 
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It is important that in the considered interval of incidence angle, the obtained phase diagram in 
figure [3] is topologically similar to the experimental one. However, nonlinear KS equation (pQ) with 
parameters defined by equations (H])-® does not presume a stable smooth surface because v y is a 
negative quantity. 




Figure 4. Relative number of islands at e — 2 keV and characteristic snapshots at = 50° (top) 
and = 63° (bottom). 



To prove that holes and ripples are stable in time, let us consider the dynamics of the surface 
morphology change. We analyze two representative kinds of patterns shown in figure [3] as snapshots 
A and E and compute the number of islands for each pattern in time. To this end, we have cut 
the surface h(x,y) at an average height level (h) and calculated the relative number of islands 
A/yiV max at fixed times, where 7V max is a maximal value of islands. In our computation scheme we 
used the following definition for the island: all points on the surface with h < (h) belonging to 
one manifold having a closed boundary, form an island. The corresponding boundary of the island 
was obtained according to the percolation model formalism. Results for relative number of islands 
were averaged over 20 independent runs. Typical evolution of the number of islands is shown in 
figure H] at e = 2 keV for = 50° and = 63°. It is seen that the relative number of islands grows 
at small time interval that corresponds to processes of the formation of islands. At intermediate 
times, the relative number of islands decreases which means a realization of coalescence processes. 
It is important that in the process of ripple formation the coalescence regime is well pronounced 
(see empty circles). On the contrary, for the process of nanohole formation (filled circles), this such 
regime is only weakly observed. At large times one has a stationary behavior of the relative number 
of islands. Hence, processes of ripple and nanohole formation are stationary ones: at large time 
intervals the averaged number of islands does not change in time. Snapshots of the silicon surface 
morphology for = 50° and = 63° at t = 0, 40, 100 and 400 seconds are shown in figure H] in the 
top and in the bottom of the figure, respectively. 
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4.2. Wavelength dependence on the ion energy 

Next, let us study the wavelength dependence on the incident ion energy and on the angle of 
incidence. As it was shown earlier in the Bradley-Harper theory, a relation between parameters 
v x and v y determines the orientation of surface patterns. The wavelength of selected patterns 
in the corresponding direction is defined as follows: A XyV = c li\^j2Kj\v m ; 1 ^ x ^ where v m i nx ,y = 
min(^ : u y ). One needs to note that following the phase diagram shown in figure [3j a variation 
in the ion energy at fixed angles of incidence causes a change in the orientation of structures: at 
small e one has structures elongated in y direction, whereas at large e, structures are horizontally 
elongated. Corresponding dependencies of the wavelength versus ion energy at fixed values for 
incidence angle are shown in figure [5j It is seen that the wavelength decreases with the ion energy 
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Figure 5. Dependence of the wavelength on the ion energy (the scaling exponent S dependence 
on the incident angle is shown in insertion). 

growth according to a power law and varies in the interval from 100 nm to 1 /im. This result is in 
good correspondence with experimental data for sputtering of the silicon target by Ar + ions [5]. 
It is principally important that as far as the penetration depth depends on the ion energy in a 
nonlinear manner (see the insert in figure [T]) one can expect a deviation from the Bradley-Harper 
wavelength asymptote A ~ e~ x l 2 . In figure [5] it is seen that at small and large incidence angles 
(see dot and dash lines, respectively) one has linear dependencies in log- log plot characterized by 
the corresponding unique slope. However, at intermediate values for related to the dash curve in 
figure [3j the dependence A(s) has a kink. This kink means a change in the orientation of patterns. 
In such a case one has two slopes at small energies, i.e., before kink, one has selected the patterns 
characterized by A x , whereas at large energies the patterns are defined by A y (see solid line and 
asymptotics in figure [5]). Therefore, at small ion energies the patterns are oriented in y direction, 
whereas at large ion energies they are oriented in x direction. Hence, for the wavelength dependence 
on the ion energy one can write A ~ e~ 6 where the scaling exponent 5 is defined as a slope of the 
dependence A(e) in double logarithmic plot before and after the kink. The dependence of the scaling 
exponent versus incidence angle is shown as an insert in figure [5] One can see that for the described 
interval for the angle of incidence 5 > 1/2. Moreover, at small and large the exponent S does not 
essentially change it values, whereas in the interval for when A(e) has a kink, the exponent S varies 
from 0.65 toward 1.05. We should note that the obtained picture is realized when the incoming ion 
flux J and temperature T are constants. In the opposite case, variation in J and T leads to the 
known asymptotes: A oc J -1 / 2 , A ex T -1 / 2 exp(— E a /2T), where E a is an activation energy. 
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4.3. Scaling properties of patterns 

Finally, let us study the scaling properties of the surface patterns, computing growth and 
roughness exponents. To this end, we analyze a height-height correlation function Ch(r,t) = (\hjr+ 
r',t) — h(r' \t)] 2 }. In the framework of dynamic scaling hypothesis following references [21, 22], one 
arrives at scaling relations Chit) oc t 2 ^, C^(r) oc r 2a , allowing one to define the growth exponent 
j3 and the roughness exponent a. 

In reference [23| it was shown that there is a set of exponents {(3} describing the universal 
behavior of the correlation function at early stages of the system evolution. At late times where 
a true scaling regime is observed there is a unique value for f3. The roughness exponent a takes 
similar values at different time windows and can be considered as a constant depending on the 
system parameters only. From practical viewpoint, the analysis of the surface growth is urgent at 
large time intervals where the true scaling regime is observed and there is no essential difference 
in values /3 at different time windows. It is known that anisotropic surfaces studied in this paper 
may exhibit a more complex dynamic scaling behaviour than isotropic ones because anisotropy 
of the surface is reflected in lateral correlations of the surface roughness 24, 25j. In reference (6] 



it was proposed to use local roughness scales a x , a y in the directions normal and parallel to the 
projection direction of the ion beam. Values for growth and roughness exponents together with 
surface tensions v x and v y at e = 2 keV and fixed values for incidence angle are presented in table [TJ 
It is seen that when v x < u y , a relation a x > a y is realized due to orientation of the structures in 
y direction. On the contrary, if v x > v y holds, then one has a x < a y . Hence, making an analysis 
of the obtained scaling exponents, one can conclude that if structures are oriented in y direction, 
then roughness is larger in x-direction and vice versa (compare patterns in snapshots A, C, D and 
E in figure [3] with exponents in table [T]). The obtained results for growth and roughness exponents 
are in goo d correspondence with experimental studies of the silicon target sputtered by Ar + ions 
(see Um^). 



Table 1 . Growth and roughness scaling exponents at e = 2 keV. 
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5. Conclusions 

Two-level modeling for nanoscale pattern formation on silicon target induced by Ar + ion sput- 
tering has been reported. We have used Monte-Carlo simulations and a continuum approach based 
on the Bradley-Harper theory. It was shown that for the described system, the dependencies of 
the averaged penetration depth of the incident ion and the corresponding distribution widths of 
the deposited energy in directions parallel and perpendicular to the incoming beam versus ion 
energy are of the power-law form. Varying the incoming ion energy and ion incidence angle, we 
have defined the sputtering yield with the help of Monte-Carlo simulations. The obtained results 
have been used in the modified Bradley-Harper theory within the framework of two-scale modeling 
scheme. 

We have computed a phase diagram for control parameters: i.e., incidence angle and ion energy 
that defines possible patterns on silicon target sputtered by Ar + ions. It was shown that at small 
incidence angles, nanohole patterns are realized, whereas at large incidence angles, pattern type of 
ripples is observed. Analyzing the morphology change of silicon surface we have shown that during 
the system evolution, the number of nanoholes/ripples becomes constant, indicating stability of 
the obtained structures in time. 
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We have found that there are deviations from the Bradley-Harper asymptotics for the wave- 
length dependence on the ion energy. Moreover, when the orientation of patterns changes, a kink 
is realized in such asymptotics. The exponent of such power-law asymptotics depends on the angle 
of incidence. At fixed values for incidence angle one has two scaling exponents related to small 
and large values for the ion energy according to a change in the orientation of structures. While 
studying the scaling characteristics of the height-height correlation function, the growth exponent 
together with longitudinal and transverse roughness exponents are obtained for different values of 
incidence angle at a fixed ion energy. It was shown that relations between roughness exponents are 
defined through relations between corresponding effective surface tensions. 

The results obtained in a two-scale modeling scheme are in good correspondence with the known 
theoretical and experimental data for sputtering of silicon target by Ar + ions in the considered 
interval of values for incidence angle of ions, the incoming ion energy, temperature and ion flux 0, 
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Morphology change of the silicon surface 



3Mma MopcfccmoriT noBepxHi KpeMHirc> npi/i po3ni/uieHHi fioro ioinaMi/i 
aproHy 

B.O. XapneHKO, fl.O. XapneHKO 

iHCTi/uyT npi/iK/iaflHoT <4di3i/i ki/i HAH YKpaTHi/i, By/i. l~leTponaB/iiBCbKa 58, 40030 CyMi/i, yKpa'ma 

npoBOflHTbC5i TeopeTUHHe floc/iifl>KeHH5i npoi^eciB 3Mim/i MopcfDO/ioriV noBepxHi KpeMHiK) npi/i po3nn/ieHHi v\oro 
ioHaMi/i aproHy b paMKax flBopiBHeBoT cxeMi/i, lu,o BpaxoBye MeTOfli/i MoHTe-Kap/io Ta MOfli/icfDiKOBaHy TeopiK) 
Epefl/ii-Xapnepa. OTpi/iMaHO Ta npoaHa/ii30BaHO c|Da30By fliarpaiviy y n/iouj,nHi KyT na^im-i* Ha/iiTaK)Horo ioHy 
Ta eineprifl ioHy, lu,o i/itocrpye MO>K/ii/iBi Ti/ini/i noBepxHeBi/ix HaHO-crpyKTyp. OTpi/iMaHO y3ara/ibHeHy CTeneHeBy 
3a/ie>KHicTb flOB>Ki/iHi/i xBi/i/ii OTpi/iMaHi/ix noBepxHeBi/ix CTpyKTyp Bifl eHepriV Ha/iiTafOHi/ix ioHiB. ripoaHa/ii30BaHO 
noKa3Hi/iK pocTy i noB3flOB>Ki-iii7i Ta nonepeHHi/m noKa3Hi/iKi/i ixiopcTKOCTi OTpi/iMai-ii/ix noBepxoHb. 

K/iK)HOBi c/iOBa: ioHHe po3niAJieHHn, Mopdponorin noBepxHi, HaHO-crpyKrypiA 
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